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Abstract. In this short note, we give an algorithm to get an explicit expression of the 
Luroth invariant in terms of the Dixmier-Ohno invariants. We also get the explicit factorized 
expression on the locus of Ciani quartics in terms of the coefficients. Finally, we answer two 
open questions on sub-loci of singular Luroth quartics. 



1. Introduction 

A Luroth quartic is a plane quartic C containing the ten vertices of a complete pentalateral, 
which is to say that there exist 5 lines £i, . . . such that Ut^j((^» = 0) ^ (£j = 0)) consists of 
10 points contained in C. Over the complex field, the closure of the locus of Luroth quartics 
in P 14 is known to be an irreducible hypersurface given by an SL 3 (C)-invariant equation, 
called the Luroth invariant and denoted L in the sequel. The classical study of the Luroth 
invariant culminated in 1919 with the work of Morley [7], which showed that the degree of L 
is 54. More recently, after the seminal work of [1], several authors have revived the subject 
in [2, 11, 10, 9] (see also [12] on the undulation invariant). However, an explicit expression 
of L was still missing. The present note explains how we were able to compute such an 
expression. Our luck was mainly being acquainted with [8], an unfortunately unpublished 
article in which Ohno gives a complete set of generators for the invariants of ternary quar- 
tics under the action of SL3(C), completing the set of primary invariants found in [3]. An 
indirect source is however [4] and the implementation of these invariants in Magma software 
by Kohel. 

In subsequent sections, we also confirm a decomposition of Luroth invariant on Ciani quar- 
tics btained in [5, Sec. 5] and answer two questions raised in [11] on the existence of invariants 
for sub-loci of singular Luroth quartics. 
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2. An algorithm to find an expression of the Luroth invariant 

Let ^3, Iq, ig, I12, I\5, lis, I2I1 I27 be the primary invariants for ternary quartics under the 
action of SL 3 (C) found by Dixmier in [3]. Let J 9 , J 12 , J15, Jig, J21 be the secondary invariants 
found by Ohno [8]. In both cases, the index specifies the degree of the invariant. Together 
these invariants form a complete set of generators for the ring of invariants of ternary quartics 
under the action of SL 3 (C). 

The method used to work out the expression of L is straightforward: 

(1) We take the 1380 monomials X = {/3 s , I^h, • • • , /| 7 } of degree 54 that generate the 
C-vector space of invariants of degree 54. 

(2) We generate a sufficiently large finite set Q of cardinality q of random plane quartics 
with rational coefficients. 

(3) We generate a sufficiently large finite set C of cardinality / of random Luroth quartics 
of the form 

where £\ = x, £2 = y, £3 = z,£^ = x + y + z, £5 is a line with random rational 
coefficients and q are rational coefficients. 

(4) We compute the matrix Mi = ((I(q))i e x,qeQi evaluating the monomials in X at the 
generic quartics in Q. 

(5) We compute the matrix M 2 = ((I(q))i e x,qeC evaluating the monomials in X at the 
generic Luroth quartics in C 

(6) We compute the 215-dimensional kernel Ni of M]_. This gives a basis of the homoge- 
neous relations of degree 54 that are satisfied by the invariants of all ternary quartics. 
Note that 1380 — 215 = 1165 is indeed the dimension of the vector space of invariants 
of degree 54, as can be computed from the Poincare series given in [13, p. 1045]. 

(7) We compute the 216-dimensional kernel N2 of M 2 . This gives a basis of the homoge- 
neous relations of degree 54 that are satisfied by all Luroth quartics. 

(8) A non-zero element in the complement of N2 in N\ is an expression for L in terms of 
the Dixmier-Ohno invariants. 

All these computations were done with Magma software. Over finite fields F p with a prime 
p = 2017, 10007, 100003 or even 1000003, computations can be done in less than a minute. 
However, getting the result over the rationals is more challenging. The main concern is to 
deal with matrices Mi and M 2 whose coefficients are as small as possible. So, at Step (2) 
of the algorithm, we generate plane quartics with random integer coefficients only equal to 
— 1, 0, or 1. Similarly, we restrict Step (3) to Luroth forms defined by integer coefficients 
Ci bounded in absolute value by 4. We can estimate the feasibility of the final computation 
through the Hadamard bounds for our matrices Mi and M 2 ; the quartics under consideration 
yield bounds slightly smaller than 2 200000 for Mi and 2 350000 for M 2 . As a sanity check before 
running the code over the rationals, we verify that these restrictions still allow us to obtain 
a valid result modulo small primes. Most of the time is spent at Step (6) and Step (7) of 
the algorithm, precisely 5 and 9 hours on our laptop (based on a Intel Core l7 M620 
2.67GHz processor). 

The program to get the result is available at 

http : // iml .univ-mrs . f r/~ritzenth/programme/luroth/luroth.m 
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It uses the implementation of the Dixmier-Ohno invariants in Magma by Kohel available at 

http : // echidna . maths . usyd . edu . au/kohel/alg/ index . html 

The 1.4Mb result is available at 

http : // iml .univ-mrs . f r/~ritzenth/programme/luroth/LurothInvF.m 

It is given by 1164 monomials with rational coefficients, the largest of which is a quotient of 
a 680-digit integer by a coprime 671-digit integer. Modulo 1000003, it starts as 

i 3 18 + 469313-Zfif + 710780/| + 969230if if I 9 + 374233i 3 i 6 7 i 9 + 276144lf if if 
+ 602674if if + 527614if if if + 538637i 3 if if + 392526i| I 6 if + 645841/f if if 
+ 914224/f if + 207808i 3 3 if + 31577i 3 i 6 if + 635768if + 668878i 3 15 J 9 
+ 507293if if J 9 + 318476i 3 / 6 7 J 9 + 59775if if i 9 J 9 + 581086if i 9 J 9 + 830307if if if J 9 
+ 804817/ 3 /f if J 9 + 6418if if J 9 + 578316/ 4 i 6 /f J 9 + 741618if if if J 9 + 452974if if J 9 
+ 36214if if J 9 + 522408i 3 i 6 if J 9 + 253043if J 9 + 469299if if jf + . . . 

3. LUROTH INVARIANTS FOR ClANI QUARTICS 
We call a Ciani quartic a plane quartic of the form 

4,1,22, 2 2 , j 4 , 22,^4 

ax + ox y + cx z + ay + ey z + jz . 

In [5, Sec. 5], using different techniques, Hauenstein and Sottile obtained the factorization on 
Ciani quartics of the Liiroth invariant as 

G A H 2 J 

with G,H,J G C[a, b, c, d, e, f] homogeneous of respective degree 6,9 and 12. Using our 
expression, it is easy to confirm their decomposition 

G = a- d- f ■ (adf - l/Aae 2 - 1 /4b 2 f - 1 /Abce - 1 /4c 2 d), 
H = (adf - l/Aae 2 

-l/4b 2 f + l/Abce+3/Ac 2 d)-(adf-l/Aae 2 +3/Ab 2 f + l/Abce-l/Ac 2 d)-(adf+3/Aae 2 -l/Ab 2 f+l/Abce-l/4c 2 d), 

j = _ l/49a 4 d 3 e 2 / 3 + 51/19208a 4 d 2 e 4 / 2 

- l/38416a 4 (ie 6 / + l/614656a 4 e 8 - l/49a 3 6 2 d 3 / 4 - 205/9604a 3 b 2 d 2 e 2 f - 3/38416a 3 6 2 de 4 / 2 + l/153664a 3 b 2 e 6 f+ 

The product G A H 2 J has 1695 monomials. Note that the total amount of weighted monomials 
in a, b, c, d, e and / in a degree 54 invariant is 3439. 

4. Singular Luroth quartics 

Let V be a three-dimensional vector space over C, and let P = ¥(S 4 V*) be the space of 
quartic forms on V. The vanishing locus of the Luroth invariant L defined above determines 
a hypersurface £ in the 14-dimensional projective space P that is independent of the choice 
of basis of V. We are now interested in studying the locus £ of singular Luroth quartics, 
or in other words the intersection £ n T>, where T> is the hypersurface in P defined by the 
vanishing of the discriminant invariant 7 2 7- Work by Le Potier and Tikhomirov ([6]) shows 
that 

£nv = A U£ 2 , 
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where £1 and £2 are irreducible subschemes of P of codimension 2 whose respective degrees 
as subschemes of £ equal 24 and 30. Moreover, while £1 is reduced, the reduced subscheme 
{^2)red of £2 is of degree 15. 

In [11], Ottaviani and Sernesi showed that no new degree 15 invariant vanishes on (£2)red, 
which implies that this scheme is not a principal hypersurface in £. By generating singular 
quartics in L\ and £2 and proceeding as above, we were able to verify that neither is there 
any new (apart from ^3 • l2i) degree 30 invariant vanishing on £2, which therefore is not a 
principal hypersurface in £ either. A similar conclusion holds of £1 as there is no degree 24 
invariant vanishing on L\. The Magma program to check this is available at 

http : // iml .univ-mrs . f r/~ritzenth/programme/luroth/SingularLurothInv.m 

A word is needed on how to generate singular quartics corresponding to points of £1 and 
£2. For £ 2 , we can proceed by using Remark 3.3 in [11]: we now choose the lines in Step 
(3) above such that three of them have a common point of intersection. 

For £1, things are slightly more involved, but basically it is a compilation of other con- 
structions and results in [11]. The program is available at 

http : // iml .univ-mrs . f r/~ritzenth/programme/luroth/GenerateLl .m 

We proceed as follows. 

(1) We start with 6 rational points pi, . . . ,p e in the plane that are sufficiently general 
in the sense that the complete linear system C of cubics passing through them has 
dimension 4. Choosing a basis c±, . . . C4 of C, we construct the Clebsch rational map 
c = (ci, C2, C3, C4) : P 2 — > P 3 and determine its image S, which is the vanishing locus 
of a quaternary cubic form F. Note that S will have 27 rational lines. The rational 
map c restricts to a birational map between P 2 and S, and we can use it to construct 
two skew lines /, m on 5*, for example by taking I to be the image of the line between 
Pi and P2 and taking m to be the image of the line between p\ and p^. 

(2) Next we determine one of the involutory points Q associated with the pair l,m of 
lines on S by the construction on page 1759 of [11]. 

(3) By Proposition 3.1(i) of [11], we obtain a quartic in £1 by constructing the ramifica- 
tion locus of the projection S — > F(TqS) from the point Q E S. 

All of these steps can be implemented using very little beyond linear algebra: 

(1) This part is straightforward. The final equation F can be determined by constructing 
the degree 3 homogeneous expressions in the q and applying some linear algebra. 

(2) We choose coordinates on / by taking two points l\, Z 2 E I and sending (x : y) E P 2 to 
xl\ + yl2, and similarly on m by choosing 7711,7712 E m. To determine the morphism 
/ : A —7- B explicitly in these coordinates, we choose two equations M\ = M 2 = 
defining m. Given [x : y) E P 2 , the point f(p) = T p S fl B corresponds to the vector 
space that is the kernel of the matrix whose rows are given by Mi, M2 and the partial 
derivatives of F. A generating vector for this space will be combination of mi and m-i 
with homogeneous quadratic coefficients fi(x,y), f 2 (x,y) in (x,y). The morphism / 
now corresponds to the map P 1 — > P 1 given by (/1, f 2 ). Similarly, one determines g. 

To construct the /' in [11] over the ground field, first compute the branch divisor 
D of g by deriving and pushing forward. If D consists of two rational points (x\ : y\) 
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and {x 2 : 2/2)5 then one can take 

f'(x : y) = {{y x x - Xxy) 2 : (y 2 x - x 2 y) 2 ). 

The case where the points of D are defined over a proper quadratic extension is most 
easily illustrated by using the affine coordinate t = x/y. Suppose that using this 
coordinate, D = [d] + [d] consists of two conjugate points, and suppose moreover 
that these points do not sum to 0. The /' = ((£ — d)/{t — d)) 2 of above now has 
the property that /' = 1/ '/'. But then one verifies that (df + d)/(df + d) is a 
fractional linear transformation of /' that is stable under conjugation and hence 
defines a morphism over the ground field. If d and d sum to 0, then one can use 
((1 + d)f + (1 + d)) /((l + d)f + (l + d)) instead. More invariantly, if D is defined 
by an equation of the form rx 2 + sxy + ty 2 = 0, then 

• if s = 0, then one can take 

- , \ / 2 2 2 2 \ 

f [x :y) — {rx — 2txy — ty : rx + 2txy — ty ); 

• if s 7^ 0, then one can take 

f'(x : y) = [r 2 sx 2 + 2r(s 2 — 2rt)xy + (s 3 — 3rst) : r{rsx 2 + Artxy + sty 2 )). 

Finally, we have to show how, given the morphisms /, /' : P 1 — > P 1 , we can compute 
a q G P 1 such that the fiber of / over q is also a fiber of /' over a point q' . For this, 
write 

/(sc : y) = (a x x 2 + b x xy + Ciy 2 : a 2 x 2 + b 2 xy + c 2 y 2 ) 

and 

f'{x : j/) = (aix 2 + b[xy + qj/ 2 : a' 2 x 2 + 6' 2 a:y + c 2 j/ 2 ). 

If we let q = (Ai, A 2 ) and = (A^, A' 2 ), then this question reduces to finding solutions 
of the equation 

A 2 (aia; 2 + b x xy + c\y 2 ) - Xi(a 2 x 2 + b 2 xy + c 2 y 2 ) 
= X' 2 (a[x 2 + b[xy + c[y 2 ) - \[(a 2 x 2 + b' 2 xy + c 2 y 2 ), 

which corresponds to the determination of the kernel of the matrix 

(ai —a 2 —a[ a 2 \ 
61 -62 -6i 6' 2 • 
ci -c 2 -ci c' 2 J 

(3) We suppose Q = (1 : : : 0) for simplicity. We get an induced coordinatization 
of TqS by sending (x : y : z) to the tangent direction given by the line through 
(1:0:0:0) and (1 : x : y : z). Then a point (x : y : z) is a ramification point of the 
projection 5 — )■ P(Tq5) if and only if the equation F(l,xt,yt, zt) = has a double 
root outside 0, or in other words if the discriminant of the quadratic polynomial 
F(l,xt,yt, zt)/t vanishes. This discriminant is a homogeneous quartic form in the 
variables x, y, z, which defines the quartic in L\ that we were looking for. 

Having generated a sufficiently large database 1 of curves in L\ by choosing random 6-tuples 
{pi, . . . ,Pe}, we can again proceed as in Section 2, which quickly yields our results. 



10,000 curves over Q are available at http : //iml .univ-mrs . f r/~ritzenth/programme/lurotri/LlDat abase .m 
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Remark 4.1. We also tried to generate quartics in C\ by using Remark 3.4 of [11]. We take 
cubics S of the form 

t 2 x + t(ax 2 + 2bxy + 2cxz + dy 2 + 2eyz + fz 2 ) + g(x, y, z) 

with g a random degree 3 homogeneous polynomial, such that S is non singular and e 2 = df . 
The last condition ensures that p = (0 : : : 1) belongs to the Hessian H of S and, after 
checking that p is non singular, we take quartics which are tangent plane sections of H at 
p. Unfortunately, it seems that these quartics are special in Li since there are degree 24 
relations between their invariants (there is a 27 dimensional space of relations in degree 24 
between randomly generated quartics of this form). 

5. Open questions 

The expression L of Liiroth invariant that we found depends on several arbitrary choices 
that may explain its cumbersomeness. First, there is the choice of the basis of invariants. 
Though some of the Dixmier invariants have geometrical interpretations that are 'natural', 
the same is far from evident for the new Ohno invariants. Secondly, our choice can be 
modified by any element of the kernel N\. Beyond a cancellation of the coefficients of 
215 of these monomials that we have already accomplished by simple linear algebra, further 
minimization of the number of monomials in the expression for L could in theory be achieved 
by techniques based on coding theory. Still, the parameters seem too large to make this 
feasible in practice. 

The negative answers on the existence of degree 24 and 30 invariants in Sec. 4 leave us clueless 
on a possible decomposition of L as the one hoped in [11, p. 1764]. 

An expression in terms of the 15 coefficients of the generic quartic would of course be 
canonical. However, it is not even practically achievable to formally express all Dixmier- 
Ohno invariants, since these contain too many monomials, as is for instance the case for 
the discriminant. A count of weighted monomials in 15 variables for degree 54 invariants 
leads to a total of 62 422 531 333. Of course only a fraction of these monomials may occur 
in the final expression of L, but we could not figure out their number, let alone the Newton 
polytope of L. 
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